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In this work a short overview of the development of spin glass theories, mainly long and short 
H ' range Ising models, are presented. 

1. Prologue 

It has been a long and hard way to unravel the fascinating subtleties involved in the physics of 
spin glasses since the pioneering experimental work of Cannella and Mydosh H in the dilute metallic 
alloy CuMn with 0.9% Mn. Nowadays these systems comprise a large variety of distinct materials 
emboding the two basic ingredients: frozen disorder and antagonistic interactions (frustration). For 
a review see for instance Rammal and Souletie Q, Binder and Young ||, Chowdury and Mookerjee 
{§, Mezard et al f|, Bray {§ and Fisher and Hertz §. 

The most successful and popular theoretical model to describe the physical properties of spin 
glasses is the one introduced by Edwards and Anderson S whose mean field version was proposed by 
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Sherrington and Kirkpatrick (SK) U . Even today some aspects of the SK model remains elusive such 
as the structure of the free energy barriers J10| , the ordering field of the condensed phase jll], [12| and 
its dynamical properties |13| . 

In this work a short account of the theoretical development in spin glasses mainly the Ising model, 
is presented. There is no intention of completeness in this work nor to give any details of the model 
calculations or a complete list of references but only of providing useful informations concerning some 
results up to now. In section 2 the mean field SK model and the picture arising from its solution is 
discussed. In Section 3 the counterpart of some special short range models where exact solutions were 
obtained are considered as well as some scaling and renormalization group theories. Some concluding 
remarks are presented. 

2. Mean field theory: the unfolding of complexity 

The much referenced mean field theory of fcrromagnetism due to Weiss may be obtained through an 
exactly solvable model. It consists of a model where all spins interact among themselves with vanishing 
size dependent interaction p4, n5| . Its solution reveals that the mathematical mechanism responsible 
for the phase transition occuring in the model is the same as in the more palatable two-dimensional 
Ising model, i.e., asymptotic degeneracy of the largest eigenvalue of the transfer matrix associated 
with the partition function of the system. It is then possible to formulate the mean field theory of 
ferromagnetism within an aesthetically atractive way as the solution of a long-range ferromagnetic 
model. 

In this same spirit, a long-range model intended to represent the mean field theory of spin glasses 
was introduced by Sherrington and Kirkpatrick || on the footsteps of the wide general Edwards and 
Anderson model ||. It is defined by the Hamiltonian 

U = - ^^ JijViVj - H X! a% W 

{id) ' 

where cr, = ±1, i = 1, 2, ...N is a set of Ising variables under an external magnetic field H. The set 
of exchange interactions coupling constants {Jij} are independent random variables chosen from the 
gaussian distribution 

p < j «)=\/S^^(-#) (2) 

and the sum is taken over all pairs (i,j) of spins. The scaling of the variance as J 2 /N is necessary 
in order to have a finite free energy per particle as N — > oo. For a system with frozen-in (quenched) 
disorder the free energy is given by 

f{Jii} = ^ L = - } §ZnZ{J ij } (3) 

for a given set of {Jij}. For a very large system it is expected that f{Jij} and others densities of 
extensive quantities are sample independent, i.e., they are self-averaging. Thus, instead of calculating 
/ for a given set (sample) as in (3), one may obtain an average free energy where the random variables 
are eliminated by carrying out the averaging of (3) over it distribution, namely 

kT 

f = -—<£nZ{J lJ }> J (4) 

where < ■■ ■ >j means this average procedure. 

It is reasonable to expect the model (1) to have many ground-states, with a complex phase space, 
in addition to the presence of many metastable states. For such a disordered model it seems almost 
miraculous that below its critical temperature T c — J/K the complex structure and organization of 
the phase space could have its details worked out. 



Historically, two complementary approaches were undertaken in order to calculate the free energy 



F{Jij} 






r 



of the model. The first was to work out directly eq.(4) through the replica method J8|, 
second to obtain F{Jij} in eq.(3) in terms of the local mean magnetizations of the spins [n6 
TAP]. 

In the TAP approach F{Jij} is given by 



, and the 
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- J^[(l + m t )£n(l + rm)/2 + (1 - m t )£n(l - m l )/2] - J^ Hm * 
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where rm =< &i > is the mean spin on the i — th site, given by 
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There are many degenerate solutions to equation (6) with the same free energy density yielding the 
ground states of the system in addition to a huge number of metastable states. While the degeneracy 
of the former solutions do not contribute to the entropy density the latter behaves like exp[Nuj(f)] 
where w(/) > for / larger than a given critical free energy density f c pj], n8l n9l |j. One has then 
the picture that below the critical temperature the phase space of the SK model has many distinct 
thermodynamic phases (or pure states) separated by infinite free energy barriers. Unlike a ferromagnet 
whose twofold degenerate ground states are related by time reversal symmetry, there is no obvious 
symmetry among the pure states of the SK model. Each state may be characterized by the set of 
local average spin < er^ > e = mf, denoting a possible equilibrium thermodynamic state. Numerical 
solution of the TAP equations |2Q] suggests that below the critical temperature, as the temperature 
decreases there occurs a continuous bifurcation (or better, a multifurcation) cascade in the number 
of solutions, a picture suggested earlier for the condensed phase reflecting the critical character of the 
spin glass phase [gl], £2], ||, || . 

The distinct thermodynamic states in which the phase space may be decomposed will have a free 
energy density f e given by g| 



exp (-Nf e /kT) = Z e = J2 exp{-E x /kT) 



(7) 
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where the sum is over all microscopic states associated with state I and Zn is the partition function 
of this state. The Boltzmann-Gibbs partition function involving all states is 



Z = J2 z t = J2 e M-Nfi/kT) 



(8) 



and a given thermodynamic state 



has a statistical weight 
Pi = \ exp(~Nf e /kT) 



(9) 



Although the free energy density is self- averaging as N — ► oo, different states may have distinct 
weights due to fluctuations in / which are 0(1/N) yielding distinct values to (8). The Boltzmann- 
Gibbs average of an observable A may thus be written 



< 



A > T = ^Pe < A> 
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(10) 



where < A > T ' is the thermal average of A in state I. There are quantities like energy and mag- 
netization which are independent of the state (reproducible) and self-averaging while others like the 



susceptibility is not. Up to now there is no known analytical way of computing the thermal average 
of an observable in a pure state. This demands the knowledge of how to project out this state say 
through an ordering field ||llj| . Although, it may be shown the certain quantities are both sample 
independent (self-averaging) and state independent (reproducible) j|, || . 

Another approach which has been successful in working out the properties of (1) has been the so 
called replica method [§L 9|. One uses the identity In Z— lim (Z n — \)/n in equation (4), interchanges 

n— *0 

the limits N — > oo and n — > 0, considers n an integer to work out < Z n >,j and at the end of the 
calculation takes the limit n — > 0. This procedure yields the following expression for the averaged free 
energy density / 
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where the parameters q a p, ot,/3 
df/dq a p = 0, which give 



(11) 



1,2, ...n are to be determined variationally from the conditions 



<7a/3 = < CT Q cr/3 > : 



Tr{a a a/3exp(/3H n )} 



with 



Tr exp(fiH. n ) 
H n = 13 J 2 ^2 qapo a V[3 + H^a a 
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where the traces in (11) and (12) are taken over n replicas at a single site. In their original solution 
Sherrington and Kirkpatrick H considered only the solution with all q a p = q, i.e., a single, order 
parameter invariant under permutation of replica labels. This solution, however, gives a negative 
entropy at low temperatures a wrong result for an Ising model |9[ . The study of the fluctuations of 
(11) around the replica symmetric solution q a p — 1 [?, |24| revealed that this is an unstable solution 
below the critical temperature, this instability persisting even in the presence of a magnetic field or 
when the interactions have a ferromagnetic component. Thus the correct solution for T < T c must 
have broken replica permutation symmetry. It took some ingenuity to find out how to break the 
symmetry among the q a p. By generalizing Blandin |25 work, Parisi |26, |27| was able to exhibit the 
correct ansatz to solve the model. It amounted to introduce an infinite number of order parameters 
q a f3 which in the limit n — > reduces to an order parameter function q(x), xe[0, 1] with the solution 
being marginally stable throught the condensed phase |26 



29, p_7|. This marginal character seems to 



reflect the critical aspect of the condensed phase and the continuous multifurcation as the temperature 
is lowered (see Binder and Young [[| for other possible explanations of these zero modes). Within 
Parisi's ansatz the free energy density takes the following functional form 



0f[q(x)) = - 



where 



with the boundary condition 
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G{l,y) = ln[2 cosh{(3y)] 



(14) 

(15) 
(16) 



During some time the physical interpretation of q a p and the broken replica symmetry solution 
leading to an order parameter function remained misterious. Parisi | p0| showed these to have a clear 
physical interpretation in terms of the overlaps of the local magnetizations and of the probability 
distribution P(q) of these overlaps between states: 



N ^— i % 
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P(q) = ^2P a PfiS(q- qa p) (17) 

a, (3 

where now a, (3 label the possible pure states with weight P a and Pp. It can be shown that P{q) 
is not a self-averaging quantity but its average over all realizations of { Jij}, P(q), is related to q(x) 
by P(q) = dx/dq. So the inverse function x = x{q) gives the cumulative probability for having an 
overlap p. Moreover, by considering any three pure states a±, as, 0*3 and the probability P(qi, q2, 93) 
for them to have overlaps q\ = q a2 a 3 , <Z2 = Qa 3 an 93 = Qa 1 a 2 it can be shown |U| that the space of 
the pure states of the SK model is organized in an ultrametric fashion: give any three states, at least 
two pairs will have the same overlap which will be less than or equal to the third pair (for instance, 

<7i =<72 < qz)- 

It is rather dificult to work directy with equations (14)-(17). For T < T c one may resort to series 



expansion of the functional equation (11) in terms of the order parameters 26, E7l E8[. However, by 
introducing a Lagrange multiplier function P(x,y), related to the local internal field acting on the 
spins, a new free energy functional may be introduced allowing solution of the SK model for all values 
of T and H [?, |33|, |3J, ^, |3(|. It is worth pointing out that in his work Temesvari j$5| argues that 
the TAP equations and Parisi's theory may not be equivalent. Although it seems that the solution of 
the SK model has been fairly worked out some points as the free energy barrier s |1(| a nd the ordering 
field [|ll[ may deserve further investigation as well as its dynamical properties |3j 
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3. Finite dimensional systems: competing unfinished theories and bizarre lattices 

Up to now there is no generally accepted theory to describe the properties of finite dimensional 
spin glasses. 

For a uniform system like a ferromagnet, its phenomenology is relatively easy to guess and a 
detailed calculation of the critical properties may be accomplished through Wilson's renormalization 
group framework. Well above T c in the paramagnetic phase the long distant spins are decorrelated 
while below T c there occurs long-range correlation among them, typical fluctuations involving clusters 
{droplets) of correlated spins. However, even for a ferromagnet model the concept of a cluster is 
somewhat vague p8[ . 

In the spin glass case, despite the heroic effort of many people, it seems that much remain to 
be done. A sound theory a la Wilson's renormalization group [?] does not exist. The main existing 
theories are the Sherrington-Kirkpatrick model and its complex many states ultrametric space with an 
Almcida-Thouless transition line in a field for one hand and the domain wall (droplets) phenomcno- 
logical approach on the other side, which in its present form does not yield the same rich structure of 
the SK model, mainly the many states and transition in a field picture | )40|, pf . On the experimental 
side there is room to fit the findings favouring one or the other framework |ll|, |42| , [l3[ the same being 
true in Monte Carlo simulation in small samples J44J, |4^, [l6|, [|7). Nevertheless, even the very existence 
of a phase transition at finite temperature in the 3D Ising short-range case is far from being settled 

s 

However it is worth to mention that certain efforts has been devoted to investigate exactly solvable 
SG models with short range interactions in bizarre lattices in attempt to understand some aspects 
of the problem not present in the infinite-range models like the correlation length, the sensibility of 



the boundary conditions and finite size effects. One of them, the Bethe lattice, has a finite number 
of nearest neighbours and so might be expected to be closer in nature to real spin glasses than the 
SK model. However due to its thin and local treelike structure the Bethe lattice contains itself some 
pathologies: there is no loops and therefore just one path linking any pair of sites, a characteristic of 
linear systems, and a finite surface to bulk sites ratio in the thermodynamic limit. These leads to very 
subtle and sensitive properties accordingly to the chosen boundary conditions. For a full discussion 
of these points see Chayes et al j|9| and Carlson et al Q, |3l], Q ■ 

Many earlier works were done on random systems on the Bethe lattice by many authors, specially 
in Japan, that derived recursions relations to find the distribution function of the effective field in the 
± J Ising SG (for a review of these works see Katsura |p3f ). However the first mean field study of spin 
glasses on the Bethe lattice as an alternative approach to the SK model was carried out by TAP Jl6f 
where the lowest order 1/z expansion was taken on the Bethe cluster. After that Bowman and Levin 
p4fl discussed the entropy and obtain the solution in the absence of a magnetic field while Thouless |55[ 
examined this model for small magnetic field in the neighborhood of the critical point. By analizing 
the correlation between two replicas he found that a replica-symmetry-breaking transition occuring on 
the same critical curve in the HT plane as obtained by de Almeida and Thouless [?] for the infinite- 
range model. Thouless [651 notice that on either sides of the critical curve the correlation functions fall 
off exponentially with distance but with one correlation length diverging on the curve. He also notice 
that the thermodynamic averages (internal energy, magnetization and the nonlinear susceptibility) 
are smooth on the transition curve except in zero field, a distinct behavior of the cusps found in the 



infinite-range model. A formal replica method was considered by Mottishaw |56| to study this model 
showing to be necessary to break the replica symmetry on the Bethe lattice just below T c to have a 
stable solution. This suggests the existence of many coexisting thermodynamics states as occuring in 
the SK model. However this conclusion was in contrast with the Thouless one |55| who found a replica 
symmetric stable solution for zero field. The controversy was later elucidated by Lai and Goldschmidt 
|57J pointing out the role of the boundary conditions in the Monte Carlo simulation of this model. 
They found that the Mottishaw's solution holds for the case of closed boundary conditions while the 
Thouless one is valid for the open uncorrelated case. This latter case has been extensively studied 
by [|50|, [5lL |5^| f° r the case of the ± J bimodal distribution. Very recently Goldschmidt |38| using the 
cavity method (see Mezard et al [g, |39|) obtained equations for the two real replicas that includes 
an extra parameter m which describes the exponential distributions of free energies of the distinct 
thermodynamic states. These equations are more general than those of Thouless which are recovered 
in m — > limit. 

Another line of approach to study SG short range models was developed after the study of the 
spin-glass behaviour in three dimension carried on by Southern and Young p0| who succeed to show 
by using the very simple scheme known as Migdal-Kadanoff (MK) aproximation |6l|, £2| that there 
was no transition for d=2 while it occurs for d=3. Recently this approach has its interest renewed to 
study chaos exponents in SG |63L p5| . In this case chaos means that the effective coupling between 
two given spins at a distance L undergoes multiple and chaotic sign changes with the temperature. It 
is found by Ney-Nifie and Hilhorst pj, p5fl that the scaling theory for symmetric SG is characterized 
by four independent exponents, the thermal ones (y and y c ) at T = and the corresponding chaos 
exponents (£ and Cc) at the criticality. The £c (> Vc) exponent appears in the scaling laws that 
describe the chaotic temperature dependence of the renormalized couplings around the critical region. 
The condition £ c > y c is found to be fulfilled in an interval of dimensions (d c = 2.46 , d+ = 3.4) in 
the MK approximation therefore including d=3. These results are based on the calculation of the 
so called autocorrelation function which measures the sensitivity of the relative deviations from the 
mean of the renormalized couplings to small changes in the initial coupling distributions. Although 
these calculations were carried on within the Gaussian projection approximation, that is, after each 
RG step the distribution of the renormalized couplings is replaced by a Gaussian of the same mean 
and width, the results were confirmed by numerical estimates obtained by maintaining numerically 
the renormalized distributions [p5|| . 

Following another line Coutinho et al J66| study within the MK scheme the structure of local 



EA order parameter of the SG Ising model instead of look to the distributions of the renormalized 
couplings. They found an exact recursion relation for the local magnetization of the model defined 
on the diamond hierarchical lattice. These lattices (hereafter DHL) are just the lattices where the 
MK aproximation is exact for the pure Ising model. They found that around the critical temperature 
a measure constructed with the normalized local EA order parameter for a n-level DHL is a fractal 
measure. The i(a) function that characterizes how the singularities of the measure are distributed 
was numerically obtained and compared with one for the pure case pCfl . 

The i(pe) function for the SG is non-trivial around the critical point while the one for the pure 
case is non trivial only at the criticality Bq] . Furthermore the former extends over a range of values 
of the a exponent much larger than the one for the latter. This fact suggests that the structure of the 
local order parameter inside the condensed phase is much more complex and requires an infinite set 
of exponents to be properly described. On the other hand dynamical simulations for 3d Ising SG in 
simple cubic lattices carried on by Bernardi and Campbell [ p7| shows that the dynamic exponent in 
the power-law relaxation of the autocorrelation function at the ordering temperature is very sensitive 
with the distribution chosen supporting that the standard universality rules do not seem to hold for 
these systems. 

The present status of the studies in spin glasses, reflects how hard the problem is and the far 
reaching consequences of its eventual comprehension. 
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